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Abstract 

Network flow interdiction analysis studies by how much the value of a maximum flow in a network 
can be diminished by removing components of the network constrained to some budget. Although this 
problem is strongly NP-complete on general networks, pseudo-polynomial algorithms were found for 
planar networks with a single source and a single sink and without the possibility to remove vertices. In 
this work we introduce pseudo-polynomial algorithms which overcome some of the restrictions of previous 
methods. We propose a planarity-preserving transformation that allows to incorporate vertex removals 
and vertex capacities in pseudo-polynomial interdiction algorithms for planar graphs. Additionally, a 
pseudo-polynomial algorithm is introduced for the problem of determining the minimal interdiction budget 
which is at least needed to make it impossible to satisfy the demand of all sink nodes, on planar networks 
with multiple sources and sinks satisfying that the sum of the supplies at the source nodes equals the 
sum of the demands at the sink nodes. Furthermore we show that the A;-densest subgraph problem on 
planar graphs can be reduced to a network flow interdiction problem on a planar graph with multiple 
sources and sinks and polynomially bounded input numbers. However it is still not known if either of 
these problems can be solved in polynomial time. 

1 Introduction 

In this paper we are interested in minimizing the maximum flow of a network by removing arcs and vertices 
constrained to some interdiction budget. This problem is known as network interdiction or network flow 
interdiction^ One can either allow or disallow partial removal of arcs (removing half of an edge corresponds to 
reduce its capacity to half of the original value) . However the techniques and results seem not to substantially 
differ on this issue. We are interested in the case without partial arc destruction. 

The problem of finding the k most vital arcs of a flow network is a special case of the network interdiction 
problem where k arcs have to be removed such that the maximum flow is reduced as much as possible. 
Network interdiction and related problems appear in various areas such as drug interdiction [19J, military 
planning [H], protecting electric power grids against terrorist attacks [T7] and hospital infection control [5]. 

The network interdiction problem was shown to be strongly NP-complete on general graphs and weakly 
NP-complete when restricted to planar graphs [15j [19]. Different algorithms for finding exact solutions 
were proposed ([8j [131 [El [H]) which are mainly based on branch and bound procedures. In [3 a pseudo- 
approximation was presented. Earlier work include [TH^. 

However when dealing with planar graphs with a single source and sink it was shown that by using planar 
duality, pseudo-polynomial algorithms for the network interdiction problem can be constructed when only 
arc removal is allowed [15 . Two of the major drawbacks of these algorithms (apart from the fact that they 
can only be applied on planar graphs) are the restrictions that only arc removals are allowed and that the 
network must have exactly one source and one sink. Vertex removal can easily be formulated as edge removal 
by a standard technique of doubling vertices and multiple sources and sink are generally handled by the 



^The term network inhibition is also used. 
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introduction of a supersource and supersink [HIS]. However, these transformations destroy planarity and 
make it impossible to profit from the currently known specialized interdiction algorithms for planar graphs. 

In this work, we are interested in the development of pseudo-polynomial algorithms for planar graphs 
which overcome some of the restrictions of previous methods. Inspired by algorithms presented in [10| , ll4 j . we 
propose a planarity-preserving transformation that allows to incorporate vertex removals and vertex capacities 
in pseudo-polynomial interdiction algorithms for planar graphs. Additionally, a pseudo-polynomial algorithm 
is introduced for the problem of determining the minimal interdiction budget needed to make it impossible 
to satisfy the demand of all sink nodes, on planar networks with multiple sources and sinks satisfying that 
the sum of the supplies at the source nodes equals the sum of the demands at the sink nodes. 

This problem is closely related to the problem of determining if a flow network is n — k secure, i.e. any 
removal of k of its components does not infect the value of the maximum flow, and can be seen as a special 
case of network interdiction. 

It is not known if network interdiction on planar networks with multiple sources and sinks is a strongly 
NP-complete problem. We show that the A:-densest subgraph problem on planar graphs, a problem for which 
neither a polynomial algorithm is known nor is it known if it is NP-complete, can be reduced to a planar 
network interdiction problem with polynomially bounded numbers as input. 

The paper is organized as follows. We begin by giving some definitions and notations in Section [2] In 
Section [3j different known complexity results on network interdiction are stated and we show how the k- 
densest subgraph problem on planar graphs can be reduced to a planar network interdiction problem with 
small input numbers. Section [4] presents a simple extension of currently known algorithms for network 
interdiction problems on undirected networks were only arc removals are allowed to the case of directed 
networks. We present in Section [5] a pseudo-polynomial algorithm for network interdiction on planar s-t flow 
networks allowing vertex interdiction. In Section [6] we present a pseudo-polynomial algorithm for a special 
network flow interdiction problem with multiple sources and sinks. 

2 Definitions and Notations 

Let (V, E) be a directed graph where V is the set of vertices, E is the set of arcs and for every arc e ^ E, 
u{e) G {0,1,2,...} denotes its capacity. Two special nodes s,t e V,s ^ t designate the source node 
respectively the sink node (the generalization to multiple sources and sinks is straightforward). We call the 
network G = {V, E,u, s,t) a flow network. For V',V" C F we denote by {V',V") the set of all arcs from 
V to V" . Furthermore, for V' C V wc denote by uj^{V') respectively u!~{V') the set of all arcs exiting V' 
respectively entering V, i.e., Lu+(y') = (V, V \ V) and uj-{V') = (F \ V) ^ For any subset V of V 
we denote by [V, V \ V] the cut defined by V . The value of the cut [V, V \ V] is J2e€u^+{V') ^(^)- I" t^e 
more general setting when every arc e € E has an additional lower bound /(e) on the capacity, the value of 
the cut [V, V \ V] is defined by v{[V' , V \ V']) = Ee^^+iv) "(e) " Eee<.-(y') Ke) f\ A cut [V\ V \ V] in 
G is called elementary if the subgraph of G induced by V is connected. For two distinct vertices s,t eV, a, 
cut [V',V\ V] is called an s-t cut if s e V , t <^ V . 

A function f : E ^ R is called a flow in G (or simply flow if there is no danger of ambiguity) if it satisfies 
the following constraints: 

i) < /(e) < it(e) y eeE 

ii) Eee.+ (.)/(e)-i:ee.-(.)/(e)=0 yveV\{s,t} 

iii) T,e<^uj+{s) /(e) - EeGa.-(s) /(e) > . 

For an s-t flow / we define its value ^{f) by J2e£uj+{s) /(e) ^ 12eeuj-{s} /(^)- A maximum s-t flow is an s-t 
flow / with maximum value. The value of a maximum s-t flow in a flow network G is denoted by p^'^^^G). 

In the context of network interdiction, for every arc and node of the network p £ VU E an interdiction cost 
c{p) G {1, 2, 3, . . . } U {oo} is associated (with c(s) = c(t) — oo). The network G = {V, E, u, s, t, c) is called 
an interdiction network. An interdiction network has unit interdiction costs if c(p) G {1,cxd} V p G V U E. 



^We also use the notation uIq and uIq to specify that we are working on the graph G. 
^Sometimes the notation is used to specify the network. 
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The network interdiction problem asks to find a set R C VU E respecting a given budget constraint c{R) := 
^^g^c(r) < B (with B E {0, 1,2, . . . }), and among all these sets minimizing the value of a maximum s-t 
flow on the graph G\ R, which is the subgraph of G obtained by removing the arcs and vertices contained 
in R (when removing a vertex, all arcs adjacent to this vertex are removed too). The value of this minimum 
maximum s-t flow corresponding to budget B is denoted by ^^^^{G) (we therefore have i/™'*'^(G') = i/^°'^{G)). 
In this context a set R C VUE satisfying the budget constraint will be called an interdiction set. An optimal 
interdiction set R minimizes the maximum s-t flow with respect to the given budget. Furthermore an optimal 
interdiction set R is called minimum if its interdiction cost c(i?) is minimum among all optimal interdiction 
sets and it is called minimal when removing any arc from the interdiction set results in a non-optimal 
interdiction set. 

We deflne the network flow security problem to be the problem of finding the minimal budget necessary 
to decrease the maximum flow by at least one unit, i.e., min{_B e {0,1,2,...} | i'g^^{G) < j/™'*'^(G')}. 

The above deflnitions and problems can easily be extended to interdiction networks with multiple sources 
and sinks with fixed supply /demand. In this case an interdiction network is given by G = (V, E, u, S, T, c, d) 
where S,T (ZV with S* n T = are the set of sources respectively sinks and the function d : V Z is the 
demand/supply function and satisfies d{s) < V s G S", d{t) > V i £ T and d{v) = OW v e V \ {S U T). 

To simplify notations a circuit C in G will be represented by the set of arcs it contains. For further 
graph-theoretical terms used in this paper and not further specified in this section we refer to [TT] . 

3 Complexity 
3.1 Previous results 

We associate the following natural decision problems to the network interdiction problem respectively the 
network flow security problem. 

Problem 1 (Decision version of network interdiction problem). Given an interdiction network G, some 
interdiction budget B E {0,1,2, .. .} and a value K E {0, 1, 2, ... }, decide whether j/|^"^(G) < K? 

Problem 2 (Decision version of network flow security). Given an interdiction network G and an interdiction 
budget B E {0,1,2,... }, decide whether z^^^''(G) < i/"^'^'^(G) ? 

It is easy to observe that Problem [2] is a special case of Problem [l] by choosing K = v'^^^{G). Otherwise, 
when working on a class of interdiction networks with a single source or sink. Problem [l] can be reduced to 
Problem |2] by the following simple construction. Suppose we have a single source s (the case of a single sink 
is analogue). We introduce a new vertex s' which replaces s as source and add an unremovable arc from s' to 
s with capacity equal to K. Problem [2] on the modified interdiction network is then equivalent to Problem [T] 
on the initial interdiction network. 

The following theorem was shown in '19] by reducing a maximum clique problem to Problem [T| 

Theorem 1 ([19j). Problem^is strongly NP-complete even when the underlying interdiction network is 
restricted to unit interdiction costs. 

Furthermore there is a trivial reduction from the binary knapsack problem (see |7] for more information 
on the binary knapsack problem) to an interdiction problem on a graph with only two vertices [19 implying 
the following theorem. 

Theorem 2 ([19 ). ProWemjTJis weakly NP-complete on planar graphs even when restricted to a single source 
and sink. 

A pseudo-polynomial algorithm for network interdiction problems on planar graphs with a single source 
and sink [TS] show that that this class of problems is not strongly NP-complete. 

By the reducibility of Problem [l] to Problem [2] Theorem [l] and Theorem |2] apply also for Problem [2] when 
working on interdiction networks with a single source or a single sink. 

It is not known whether the class of interdiction problems on planar graphs with multiple sources and 
sinks is strongly NP-complete. Furthermore, the direct reduction from Problem [T] to Problem |2] is not 
possible anymore on this class of networks. We will introduce in Section [6] a pseudo-polynomial algorithm 
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for Problem [2] on the class of planar interdiction networks with multiple sources and sinks and such that the 
sum of the demands is equal to the sum of the supplies and equal to the maximum flow in the initial network. 
This algorithm does not seem to generalize in a simple way to Problem [T| 

3.2 Relation between planar network interdiction and the A;-densest subgraph 
problem in planar graphs 

We will show that finding dense subgraphs of a given size on planar graphs can easily be modelled as a 
planar network interdiction problem. The problem of finding a densest subgraphs of size k is often called the 
fc-densest subgraph problem or the fc-clustering problem and is formally defined as follows. 

Problem 3 (fc-densest subgraph problem). Given an undirected graph G = (V,E), find a subgraph over k 
vertices with a maximum number of edges. 

Whereas Problem [s] is known to be NP-complete on a wide variety of graph classes U , its complexity for 
the class of planar graphs is still open. A slight modification of Problem |3] obtained by imposing that the 
subgraph must be connected was shown to be NP-complete on planar graphs 9 . 

Theorem 3. The k-densest subgraph problem on planar graphs can be reduced in polynomial time to a 
network interdiction problem on planar graphs. 

Proof. Let G = {V, E) be a planar undirected graph. Consider the following planar interdiction network 
G' = {V , E' , u, S, T, c, d). The underlying graph {V , E') is obtained from G by subdividing all edges, i.e., on 
every edge e & E, a, new node is added. We thus obtain a bipartite planar graph where each edge has one 
endpoint in V and the other one in the set of newly added vertices Ve {V — V^Ve). By directing all edges 
from V to Ve, we get E' (c.f. Figure [T]). The sets containing the sources and sinks are defined as follows 
S ^ V , T = Ve, all arcs have unit capacity u(e) = 1 V e € £", sources have infinite supply d{s) = oo W s € S 
and every sink has unit demand d{t) = 1 G T. Furthermore all arcs and all vertices of Ve are unremovable 
(they have an interdiction cost of oo) and the vertices in V have an interdiction cost equal to one. 

For some fixed budget B e {0,1,2...}, an optimal interdiction set G' corresponds exactly to the vertices 
of a -B-densest subgraph in G because of the following observation. For some interdiction set R, the decrease 
of flow by removing the components in R corresponds to the number of sinks for which both neighbors are 
in R. This corresponds to the number of edges in G that have both endpoints in R. □ 

4 Planar duality and current pseudo-polynomial algorithms 

Planarity is a very helpful property when dealing with interdiction problems as the problem seems to have 
a simpler form when restated on the planar dual of the original interdiction network. We flrst introduce the 
planar dual of an interdiction network, which can be seen as a generalization of the classical planar dual. In 
a second step we propose a pseudo-polynomial algorithm for planar network interdiction with a single source 
and a single sink and without vertex removals. This algorithm is a direct generalization of an algorithm 
introduced in [15 , which was designed only for undirected networks]^ The extensions we propose in the 
following sections will overcome some restrictions of this algorithm. 

4.1 Planar duality for interdiction networks 

The classical planar dua|^ of a directed graph is constructed on the base of a planar embedding by placing 
a vertex in each face of the original graph and connecting two vertices by an arc if they correspond to faces 
in the original graph sharing an arc. This gives a natural one-to-one correspondence between arcs in the 
original graph and arcs in the dual graph (dual arcs) as well as faces in the original graph and vertices in the 
dual graph, and vice versa. By convention, the dual arcs are oriented such that they cross the corresponding 
original arcs from right to left. See [T^ for more details. 

* Furthermore the algorithm we present does not allow partial arc removals whereas the algorithm presented in |15| did allow it. 

This makes no real difference as the technique applies easily to both cases. 
^It is also called geometric dual or simply dual. 
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{V,E) {V',E') 
Figure 1: Topology of the auxiliary graph {V',E') used for proofing Theorem |3] 



For our purposes we extend the notion of planar duality on networks with lower and upper capacities on 
the arcs and interdiction costs on the arcs. Let G = {V, E,l,u,c) be a directed planar network where for 
every arc e £ E, l{e) , u{e) , c{e) G {0, 1, 2, ... } correspond to the lower capacity bound, upper capacity bound 
and interdiction cost of arc e (where ^(e) < u{e) Ve G i?). We define the dual G* = {V* ,E* ,\* ,c*) of the 
network G in the following way. The graph {V* ,E*) is the planar dual in the classical sense of the graph 
(V, E) with the single difference that for every arc in the dual we added a reverse arc. For every arc e € E we 
denote by the corresponding dual arc (as in the classical sense) and by its reverse arc (cf . Figure [2]) . 
The function A* : E* ^ Z is an integral length function in the network G* , defined by A*(e^) = u{e) and 
A*(eg) = -l{e) ye&E. The cost function c* is defined by c*(e^) = c(e),c*(eg) = 0^6 e E. 

For every cut [V',V \ V] in the original network we denote its corresponding dual arcs by C*{V') = 
{e'^ e E* \ e e {V, V \ V')} U {ej^ e E* \ e e {V \ V, V')}. Note that the set C*{V') is a set of edge- 
disjoint, counterclockwise non-overlapping circuits in (V*,E*), where non-overlapping is defined as follows. 
Let CI , C2 be two circuits in G* and Vi , V2 C be the vertices in V surrounded in counterclockwise sense 
by CI respectively C2 ■ We say that CI , C2 do not overlap if Fi fl V2 = 0. The following proposition highlights 
the correspondence of minimal cuts in the network G and sets of non-overlapping, edge-disjoint circuits in 
its dual G* . 

Proposition 1. The function that associates with every cut [V' ,V\V'] its corresponding dual arcs C* (V') is 
a one-to-one mapping between cuts in G and sets of non-overlapping, edge-disjoint, counterclockwise circuits 
in G* . Furthermore, the value of a cut in G is equal to the sum of the lengths of the corresponding dual arcs 
in G^ , i.e., for any subset V C V, we have 

V\V'])= ^*(^*)- 
e*ec'(y') 

Proof. The one-to-one property follows easily by observing that for any set Cq of non-overlapping, edge- 
disjoint, counterclockwise circuits, the set V' of all vertices being surrounded in counterclockwise sense by 
one of the circuits in Cq satisfies C*{V') = Cq. The equality between the value of a cut in G and the sum of 
the lengths of the corresponding dual arcs follows directly from the definition of A* . □ 

Proposition [1] implies that for every minimal cut in G, there exists a corresponding circuit in the dual G* 
with length equal to the value of the cut. In particular, when working on a fiow network G = (V, E, I, u, s, t) 
with a single source s and a single sink t, every minimal s-t cut in G corresponds to a counterclockwise 
circuit around s with t separated from s by the circuit and with length equal to the value of the cut. We call 
circuits that separate s and t in such a way s-t separating counterclockwise circuits. In the following we see 
how this correspondence can be exploited for solving network interdiction problems on planar graphs with a 
single source and a single sink. 
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Figure 2: Example dual graph {V*,E*) drawn over the given original graph 
iV,E). 

4.2 A pseudo-polynomial algorithm for single source, single sink network inter- 
diction on planar graphs without vertex removal 

We now construct a pseudo-polynomial algorithm for solving the network interdiction problem on planar 
directed graphs with a single source s and a single sink t and without vertex removal, which is a direct 
generalization of an algorithm presented in |15j (that only worked for undirected graphs). This algorithm 
nicely illustrates the techniques currently used for creating pseudo-polynomial network interdiction algorithms 
on planar graphs. Given is an interdiction network G — (V^ E,u, s,t, c) with unremovable vertices, i.e., 
c(r) = oo V r € y , and interdiction budget B. For every interdiction set i? C -E we fix a minimum s-t cut in 
G \ R that we denote by [Vr. V\Vb]. We therefore have y'^^^{G \ R) = v{[Vr, V\Vr])- Ee^mu^+iVn) 
Note that an optimal, minimal interdiction set R always satisfies R C uj^{Vr), as otherwise the interdiction 
set R' = Rf] lo'^{Vr) would reduce the maximum flow by the same value as R and has lower interdiction 
cost. The reduced value of an s-t cut G (with respect to the budget B) is defined as the minimum value of G 
va G\R over all interdiction sets R and is denoted by vb{G). Note that the problem to find for some given 
s-t cut G an interdiction set R that minimizes the value of C in G \ i? is a binary knapsack problem. 

The main idea of the algorithm is to find an optimal, minimum interdiction set R by finding a corresponding 
s-t cut with minimal reduced value. This is done by translating the problem into the dual. For any set of 
edges U* C E* , we define its reduced length (with respect to B) by X*g{U*) — "cam.{}2e*£U'\X' -^*(e*) | 
X* C C^*, X]e*ex* c*(e*) < B}. Similarly for a walk W* in G* going along the arcs (e^, 63, . . . , e\), we define 
\*b{W*) = niin{^jg^^2 fc}\/^*('Si) I ^ ^ {1: 2, . . . , A:}, c*(e*) < B}. By the correspondence between 
s-t cuts in G and s-t separating counterclockwise circuits in G* as highlighted in Section |4.1| we have that 
the problem of finding an s-t cut in G with minimal reduced value is equivalent to finding an s-t separating 
counterclockwise circuit with minimal reduced length in the dual. Such circuits can be described in the 
following way. Let P be any path in the graph G from vertex s to vertex i, we define — {e^ EE* \ e E P} 
and Pj^ = {e^ EE* \ e E P}. For any set of edges U* C E* we define its parity with respect to P by 
Pp{U*) = \U* n P^\ - \U* n PJ^\. By a result of [H] we have that for every circuit C* in G* , two consecutive 
crossings on P alternate between left-right crossing and right-left ones. This implies that every circuit C* 
in G satisfies pp{C*) E {—1,0,1}. Another implication is that a circuit C* in G* has the properties to be 
counterclockwise s-t separating if and only if pp(C*) — 1. 

We therefore have to solve the following problem. 

Problem 4. 

argmin{A|^(C*) | C* circuit in G* with ppiC*) = 1} 
Consider the following relaxation of Problem |4] 
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Problem 5. 

argmin{A^(W^*) | W* dosed walk in G* with pp{W*) — 1} 

A solution to Problem |4] can be easily obtained on the base of a solution W* to Problem[5]by the following 
observation. W* can be partitioned in a disjoint union of circuits Cl,C2, . ■ ■ ,0^. Furthermore, by modularity 
of the parity function pp and the fact that pp{W*) — 1, we have 1 — pp{W*) = Ppi^i)- the parity of 

each circuit is in {—1, 0, 1} we have that there is some index i g {1, 2, . . . , fc} with pp{C*) = 1. From C* C W* 
follows that AJj(C*) < X*BiW*). By optimahty of W* for Problemjsjwe thus have X*b{C*) = X*p{W*) and by 
the fact that Problem [5] is a relaxation of Problem [4] follows that C* is an optimal solution for Problem [4j 
More generally, the above reasoning shows that minimal solutions of Problem |5] correspond to solutions of 
Problem m and vice versa. 

We finally show how Problem [s] can be solved by a dynamic programming approach that we realize as a 
sequence of shortest path problems on an auxiliary graph G' — {V, E' , A') with positive arc lengths A' which 
is defined as follows. Let P be a shortest path from s to t in G and we denote by |P| the length of the path 
which is the number of arcs used in P. V' consists of \V* \ ■ {B + 1) ■ {2\P\ + 1) vertices that we denote in 
the following way V = \v* GV*,be {0, 1, . . . G {-\Pl -\P\ + !,■•■, \P\}- The set of arcs E' is 

defined by E' = E'f^U E[ E[ E[ _ U E!^ E'^ + U E!^ _ with 

. E^ = {{vl^, v;^,^^) I ev',b& {1, 2, .'. . , B};p e + i, . . . , \P\}} 

• - {Kp,<p) I {u*,v*) e E* \ {P^ u P^), 5 e {0, 1, . . . , B},p e {-\P\, -\P\ + 1, . . . , |P|}} 

• K+ = {Kp>K,p+i) I eP^,be {0, 1, . . .,B},p e {-\P\, -\P\ + 1, . . . , |P| - 1}} 
. E[^_ = I i^*,v*) ePE,be {0, l,...,B},pe {-\P\ + i, -|P| + 2, . . . , |P|}} 

• ^2,= = {K,p,v;^p) I {u*,v*) e E*\{PDUP^)M e {c*{u\v*),c*{u*,v*) + l,...,B}M -foi- 
c*(u*,z;*),pe{~|P|,-|P| + l,...,|P|}} 

. ^^,+ = {Kp,<p+i) I eP^,6i e {c*{u\v*),c*{u\v*) + I, . . . ,B}M = bi - c*{u* ,v*),p e 

{-IPI,-IPI + 1,..'.,IPI-1}} 

• ^2,- = {Kp,«^,p-i) I {u*,v*)ePEM^{c*{u*,v*),c*{u*,v*) + \,...,B},b2^b^-c*{u*,v*),pe 
{-|P| + 1,-|P| + 2,...,|P|}}. 

We define a function -q : V' U {E' \ Eg) V U E that maps elements of G" to corresponding elements in 
G* in the following way. 

The length A' : ^ {0, 1, . . . } is defined as follows. 

fo ife' eE',UE',^^UEi+UE'2^_ 

^ ' [A*(?7(e')) li e' (=, E[_\J E[^^\J E[_ 

We now define a correspondence between walks in G' and walks in G* by extending the function rj in the 
following way. Let ul^^^^ , vl^,p^ & V , W be a walk in G' from to v^^ p^ and {e[,e'2, . . . , ej.) the suite of 

edges corresponding to the walk W' . Furthermore, let I = {i ^ {1,2, . . . , k} \ e'^ £ E' \ Eq}. The walk ri{W') 
in G' is defined to go along the edges {f]{Si))iGi- It is easy to verify that ri{W') is effectively a walk (from 
u* to V*) in G* and satisfies pp{i]{W')) — p2 — Pi- We associate to the walk W a set R{W') C E defined 
by R{W') ^ {e e E\3i € I with e- e E!^^^ U E!^^^ U E!2 _ and = r?(e-)}. By construction of G' we have 
ciRiW')) < 5i - 62 Therefore, for any path W in G', the set P(VF') is an interdiction set. Furthermore, 
we have Xl_^_ij^{r]{W')) < X'{W') because by setting the length of the arcs R{W')^ in G* to zero, the length 
of r]{W') is smaller or equal than X'{W') 

In particular, if the walk W goes from q to i for some v* S V* , we then have that ri{W') is a 
closed walk in G* with pp{t]{W')) — 1 and satisfying X*g{r]{W')) < X'{W'). Conversely, for every circuit 
C* in G* satisfying pp{C*) = 1, we can find a corresponding path P' in G' with A'(P') = A|j(G*) in the 
following way. Let (e^, ej, . . . , e|,) be a suite of arcs in E* corresponding to the circuit G* and v* G V* be 



^When ri{W') corresponds to a path or a circuit, this relation is always satisfied with equality. 
^Here too, in case that rjiW') corresponds to a path or a circuit, we have (»)(iy )) = \'{W'). 
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the vertex corresponding to the tail of e^. Let / C {1,2, ... ,k} be a set satisfying J2i<£i i^i) — ^ 

The path P' will be defined by its corresponding suite of arcs (e'j^, 62, . . . , ej./) as follows. We set k' = 
k + B — J2iei P' starts at the vertex w^.o- The path P' will be defined by adding edges step by step. 

For i e {1,2, ... ,A:} let e* = {u* ,v*) and u^^ be the endpoint of the currently constructed path. If i G /, 
we set e[ = {ul,piK-c'(ei)^p+pp(et)) ^nd otherwise e[ = {ulp,K,p+pp{e;))- For i E {k + 1, k + 2, . . . , k'}, we 
set e[ = {vl,_^^i I, i). It is easy to verify that P' is effectively a path in G" with ri{P') = C* satisfying 
A'(P') = A^(C*)'. 

Let V* G V* be some fixed vertex and P' be a shortest path from v*^ q to Vq i. The discussion above shows 
that 'i]{P') is a solution to the following problem. 

Problem 6. 

argmin{A^(VF*) | W* closed walk in G* containing vertex v* and satisfying pp(W*) = 1} 

Solving Problem [6] for all vertices v* e V* will solve Problem [5] However, as a walk solving Problem [5] will 
pass at least once by an arc e' with ri{e') G P^ it suffices to solve Problem |6] for all vertices in G* having 
at least one outgoing edge contained in P^ . Therefore, it suffices to solve \P\ instances of Problem |6] to get 
a solution to Problem [5] Notice that solving the network interdiction problem with the above algorithm for 
some budget B, solves also the interdiction problems for all budgets B' < B. This solution can be obtained 
in an analogue manner as for the case with budget B by looking in G" at walks from Vg q to v*^_b' i- 

Applying Dijkstra's algorithm for shortest paths to solve the subproblems of Problem [6] the algorithm 
introduced above to solve Problem [5] gets an overall complexity of 0{B\P\'^n\og{nB)). Because of the special 
structure of the network G" it is easy slightly improve the running time by handling the vertices in Dijkstra's 
algorithm level by level with respect to the remaining budget, i.e., we handle all vertices corresponding to a 
budget B first (always taking the one with the smallest label as usual), then all with budget B — \ and so 
on. When looking for the vertex with the smallest label, we only have 0{n) candidates instead of 0{nB). 
This allows to reduce the running time to 0{B\P\^nlog{n)). By reversing the roles of budget and length in 
the proposed algorithm, one can replace B in the above running times by iy'^^^{G). 

Notice that the proposed method can be slightly simplified (without influencing the above running time 
bounds) allowing for the parity only to take values in {— [^] , — + 1, • • ■ 1 [^1 } instead of { — |P|, — |P| + 
1, . . . , |P|} as every circuit in G* with parity equal to one contains no subpath with a parity not contained 
in this range because it would not be possible to come back to parity one without going two times through 
a same vertex. 

The formulation of Problem [6] as a shortest path problem on G" is essentially a dynamic programming 
realization of a multi-objective shortest path problem with the three objectives budget, length and parity on 
a graph G defined as follows. The graph G — {V* ,E) is obtained from the graph G* ~ {V* ,E*) by doubling 
every arc. For every arc G E* , we denote by 61,62 the two corresponding parallel arcs in E. To every 
arc in E we associate a parity value, a length and a budget value. The parity value of ei and 62 is set to 
pp{e). The length of ei is equal to u{e) and the length of 62 is set to zero. Finally, the budget value is set 
to zero for ei and to c(e) for 62. Therefore, the arcs in E indexed by one correspond to non-removed arcs 
and the ones indexed by two correspond to removed arcs. For some fixed vertex v* € V* , a, closed walk in G 
containing v, having parity equal to one, a budget value bounded by B and with minimal length among all 
those closed walks corresponds exactly to a solution to Problem [6j 

5 Incorporating vertex interdiction and vertex capacities 

Let G = (V, E,u, s,t,c) be an interdiction network. In this section we show how vertex interdiction and 
vertex capacities can be incorporated in the algorithm of the previous section. We begin by introducing the 
possibility of vertex interdiction and observe afterwards how vertex capacities can be added to the model. The 
role of the dual network G* will be replaced by a modified dual G* which allows to model vertex interdiction 
basically as arc interdiction. This technique was used in |10j for modelling vertex capacities in planar flow 
problems. As we will see in this section, this model can also be used for modelling vertex interdiction. The 
modified dual we present in this section and extend in the next one is slightly different from the one presented 
in [lOj . using some less artificially added vertices. 
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The modified dual network G* = {V* , E* ,\* ,c*) is an extended version of G* = {V* , E* , \* , c*) with 
additional vertices and arcs defined as follows. For every vertex v & V we denote by f*{v) the face of G* 
corresponding to v. Similarly, for every v* G V* we denote by f{v*) the face of G corresponding to v*. The 
network G* consists of the vertices V* and contains an additional vertex for every face of G*^Because of 
the one-to-one relation between faces in G* and vertices in G we set V* = V* U V. The arc set E* is defined 
by E* — E* U E^ where E^ contains the two arcs {v,v*) and {v*,v) for every pair oi v G V and v* G V* , 
where /(«*) is a face adjacent to v (c.f. Figm-cjs]). The length function A* is an extension of A* on the arcs 
E* giving a value of oo to all arcs in E^ (or simply a high finite value assuring that the arc is never used in 
following applications of shortest path algorithms). Furthermore, the modified interdiction cost function c* 
is an extension of c* on the arcs in E* defined as follows. 

c*{e*) = c*(e*) Ve* G E* 

c*{v,v*) =0 y{v,v*) G E^nV xV* 

c*{v*,v) =c{v) y{v*,v) G E^ nv* xV 

For some subset of arcs U* C E* we define their reduced cost in G* with respect to the budget B by 
X%{U*) := minlX^^gjj-.^x- I ^* ^ f^'Sssgjc- c*(e*) < B}. The main idea of this model lies in the 

following correspondence between interdiction sets in G and s-t separating counterclockwise circuits in G*. 

Theorem 4. 

i) For every interdiction set R d V U E with respect to the budget B, there is an s-t separating counter- 
clockwise circuit C* in G* satisfying X*g{C*) < i'"^^^(G\ R). 

ii) For every s-t separating counterclockwise circuit in G* and every budget B , there is an interdiction set 
RcVUE with respect to B satisfying iy™'^^(G \ R) < A|j(C*). 

Proof. Let i? C V U E he an interdiction set with respect to some fixed budget B and U C E he a, set of arcs 
corresponding to a minimum s — t cut in G\R. As there is no path from s to t in the graph G \ (i? U J7) , it is 
easy to see (by means of the Jordan curve theorem) that there is an s-t separating counterclockwise circuit 
C* in G* consisting only of arcs which arc cither adjacent to vertices in R or are dual arcs of arcs contained 
in R or U. We therefore have as desired A|j(G*) < i/'°^''(G \ R). 

Conversely let C* he an s-t separating counterclockwise circuit in G* satisfying X*g(C*) < oo (when the 
reduced dual length is equal to oo, the result follows trivially). Let U* C E* he a. solution of 

argmm{ ^ A*(e*) | c*{X*) < B} . 

By definition of the reduced dual length we have A^(C*) = X*{C* \U*). Let he the subset of vertices 
in V through which the circuit C* passes and let U* = U* O E* . Because of A^(C*) < oo we have that all arcs 
of C* entering in one of the vertices in are contained in U*. The cost of U* can therefore be reformulated 
as follow. 

c*{U*) = c*{U*)+c{V^.) 

Let U = {e G E \ G U*} and we define R = U U. By the above equation and the definition of U* we 
have c{R) = c*{U*) + c(V^) ~ c*{U*) < B showing that R is an interdiction set with respect to the budget 
B. Let E~ ^ C* nE* and i?^ = {e G E \ e'^ G E~}. The fact that C* is a counterclockwise, s-t separating 
circuit implies that there is no path from s to i in G \ (V^ U E^). Therefore, removing the arcs E^ \ U 
from G \ R destroys all paths from s to i implying i'™'^^{G \ R) < u{E^ \ U). The result is finally obtained 
by observing that u{E^, \U) ^ X*{C* \ (7^) = A*j(C*). □ 

Theorem |4] implies that the problem of finding by how much the value of a maximum fiow can be reduced 
through interdiction reduces to finding an s-t separating counterclockwise circuit with minimal reduced value 
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Figure 3: Topology of the auxiliary graph (V* , E*) used for modeling vertex 
interdiction. 



in G* . Furthermore, the proof of Theorem ^ shows how one can transform such a circuit to an optimal 
interdiction set. 

For characterizing s-t separating counterclockwise circuits, we introduce an adapted version of the parity 
function. As in the previous section, let P be a path in G from vertex s to vertex t. Because in the graph 
G* it is possible to cross P at a vertex, we have to take this possibility into account in the parity function. 
We therefore define a parity function pp which is an extension of pp on the subsets ojji!* in the following 
way. For every vertex v € V \ {s, t} which lies on the path P and every arc {v, v*) € E* which leaves P on 
the left side, we set pp{v, v*) — 1. Similarly, for every vertex v E V \ {s,t} on P and every arc {v* ,v) £ E* 
which enters P from the left side we set pp{v* , v) = —1. For all other edges in E* we set pp = 0. Finally, for 
any set U* C E* we define its parity by pp{U*) = Ppie*)- As in the previous section, we have, by 

this definition of pp. that a circuit C* in G* has the properties to be counterclockwise s-t separating if and 
only if pp{C*) = 1. 

Applying the techniques of the previous section to the modified dual G* allows therefore to solve network 
interdiction problems on planar graphs with a single source and sink and with the possibility to interdict arcs 
and vertices. The^symptotic worst case complexity stays the same as in the previous section because the 
size of the graph G* is only at most a constant factor larger than G* , and the same is true for the number 
of times Problem [5] has to be solved. 

Upper capacities on the vertices can be introduced in the same way as shown in |10l by slightly modifying 
the network G* as follows. Suppose that some vertex v E V has an upper capacity u(v) G {1, 2, . . . }. Let 
V* — {v* E V* I {v* , v) e E*}. For every v* € V*, an additional arc from v* to v is added with interdiction 
cost oo and capacity u(v). For a justification of this construction in the case without interdiction (respectively 
by setting P = 0), see [TU]. Theorem |4] remains true for the modified network modelling vertex capacities. 
This can be proven in an analogue way as for the case without vertex capacities. To simplify further results, 
we consider in the following only fiow networks without vertex capacities. Vertex capacities can easily be 
added by the above construction. 

The method presented in the next section to solve the network flow security problem on networks with 



multiple sources and sinks differs signiflcantly from the approach presented in Section 4.2 Therefore the 



model presented in this section does not apply directly to this case. We will see in the next section how the 
model for vertex interdiction can be adapted (the adaption for vertex capacities will be analogue). 



6 Pseudo-polynomial algorithm for network flow security 

In this section we propose a pseudo-polynomial algorithm for the network flow security problem on planar 
graphs with multiple sources and sinks where the sum of the demands is equal to the sum of the supplies 
and equal to the maximum flow. The method is inspired by an algorithm introduced in |14| for testing if all 
demand can be satisfled in a planar flow problem with multiple sources and sinks. We will therefore begin 
with a short summary of this maximum flow algorithm for planar graphs in Subsection |6.1[ Our algorithm 



for network flow security will be introduced in Subsection 6.2 restricted to the case without vertex removals 
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This restriction will be lifted in Subsection |6.3| by adapting the modelling idea of Section |5] 
6.1 Maximum flow algorithm of Miller and Naor 

In this section we briefly present an algorithm introduced in ^14] for testing if all demands can be satisfied 
in a flow problem with multiple sources and sinks such that the sum of all demands equals the sum of all 
supplies. Let G = (V, /, u, S, T, d) be a planar flow network with lower and upper limits on the capacities 
designated by I and u, with source set S <Z V and sink set T G V \ S and with demand/supply function d 
satisfying ~d{S) — d{T). A flow in G is called saturating if it satisfies all demands. The problem we want to 
solve is the following. 

Problem 7. Does there exist a saturating flow in G? 

In a first step, the problem is reduced to a circulation problem on a network G, which is the problem 
of finding a flow in a network with lower and upper capacities on the arcs but without sources and sinks. 
The network G — {V — V,E,l,u) is defined as follows. Let T' be an undirected tree over the vertices in V 
that spans the sources and sinks and that can be added to G without destroying planarity. For every edge 
{v, u} e T' we denote by Vr'{v, u) C V the set of vertices connected to vertex v mT'\ {v, u}. We orient the 
edges in T' to obtain T in the following way. For {v,u} G T' we orient the edge from w to m if d{Vr' {v,u)) > 0, 
otherwise we orient the edge from u to v. The set E is defined to be EUT. Furthermore, the lower and upper 
capacities I and u are extensions of I and u on the set E defined by u) — u{v, u) = d{Vr' {v,u)) V(i', u) G T. 
As noted in [14] we have the following theorem. 

Theorem 5. There exists a saturating flow in G <^ there exists a circulation in G. 

A circulation in ajplanar flow network can be computed by solving shortest path problem in its dual in the 
following way. Let G* = {V* , E* , X*) be the dual network of G as defined in Section [i] with the difference 
that we have no dual costs c* as we deal with a standard flow network and not an interdiction network. Let 
r* G V* he an arbitrary vertex in G* and for v* £ V* let n{v*) be the distance of a shortest path from r* to 
V* (with respect to the length A*). In [T3j the following theorem was proven. 

Theorem 6. There exists a circulation in G <^ G* has no negative circuits. 

In this case, a circulation can be obtained in the following way. Let {v,u) G E and {v* ,u*) = {v,u)^ its 
corresponding dual arc. A circulation can flnally be deflned by assigning a flow equal to max{0, fJ.(u*) —ii(v*)} 
to each arc (w, u) £ E. 



6.2 Network security on planar graphs with multiple sources and sinks 

Let G = (V, E,l,u, S,T, d, c) be a planar interdiction network satisfying —d(S) — d(T) and without vertex 
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removal. Let {V , E, I, u) be the auxiliary network corresponding to (V, E, I, u, S, T, d) as deflned in Section 
We extend this auxiliary network to G = {V , E,l,u,c) where c is an extension of c on the set E, satisfying 
c(e) = ooWe e E\E. Let G* = {V* , E* , \* ,d^) be its corresponding dual network as defined in Section |4j 
Using Theorem [g] we can formulate a problem on G* , which is equivalent to Problem [i] for G, as follows. 

Problem 8. Does there exist a circuit C* in G* with X*g{C*) < i^'"^'^(G)? 

The circuit with minimum reduced cost can be found by similar techniques as the ones presented in 
Section [4] with the differences that we do not have to take parity into account and that we have to use a 
shortest path algorithm, which can deal with negative arc lengths. A simple implementation would be first 
to determine shortest paths for all pairs of vertices in G. By using an algorithm presented in [5] this can be 
done in 0{n? \og^ n) time. To check if there is a circuit with negative reduced value in G* going through some 
fixed vertex v* G V* , a, shortest path in an auxiliary network as described in Section [Z] will be determined. 
This can be done by determining the shortest paths from some initial vertex to all others by proceeding level 
by level with respect to the remaining budget using at each level information of the preprocessing step. Over 
the B budget levels, this can be done in 0{Bn^) time. By simply performing this operation for all possible 
start vertices, we get a running time of 0{Bn^). By reversing the roles of budget and length, we can get 
an algorithm with running time 0{v'^'^^{G)n^). We expect that this running time can even be improved by 
exploiting more deeply the structure of the auxiliary graph. 
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6.3 Generalization to the case with vertex interdiction 

Let G = {V, E,l,u, S,T,d,c) be a planar interdiction network satisfying —d{S) = d{T) and allowing arc 
and vertex removal (except for sources and sinks). As in the case without vertex removal, we begin by 
reformulating the problem as an interdiction problem for circulations. Let G = {V = V, E,l^u,c) be the 



auxiliary graph as defined in Section 6.2 and as before we denote hy T = E \ E the added tree arcs. We 
now discuss how arc and vertex removal in G can be translated to G such that Theorem [s] remains valid for 



the resulting networks. As already exploited in Section 6.2 removing an arc of G corresponds to removing 
the same arc in G. However, vertex removal cannot be translated in such a direct way as the arcs in E \ E 
are auxiliary arcs which should not be removed by a vertex removal. For any interdiction set R d V U E, 
we denote by G{R) the graph obtained from G by removing all arcs contained in R and all arcs in E being 
adjacent to a vertex in R. We have the following relation. 

Theorem 7. For any interdiction set R of G we have the following equivalence. 

There is a saturating flow in G \ R there is a circulation in G{R). 

Proof. Let G{R) be the network obtained from G by removing all arcs in R and all arcs adjacent to vertices 
in R. We trivially have that there is a saturating flow in G \ R ii and only if there is a saturating flow 
in G{R). The network G{R) can easily be obtained from G{R) by applying the construction introduced in 
Section [O] Applying Theorem |5] proves finally the claim. □ 

By a classical characterization for feasibility of circulation problems (see [T) the following theorem follows. 

Theorem 8. For any interdiction set R in G we have: 

There exists no circulation in G{R) <^ there exists a cut in G{R) with value < 0. 

Furthermore, as the arcs contained in a cut can be partitioned in groups of arcs corresponding to elementary 
cuts, we have that there is a cut in G{R) with value strictly less than zero exactly when there is an elementary 
cut in G{R) with value strictly less than zero. In the following, we show how the problem of finding an 
interdiction set R and an elementary cut in G(i?) with negative value can be mapped onto a modified dual 
networ k of G. Let G* = {V* , E* , X* ,c*) be the modified dual network for the network G as introduced in 
Section sj^ Analogously as in Section [s] we use the notations V* — VLlV* and E* ~ E*\JE^ . We associate 
with every circuit C G G a modified reduced length 7|j(C) defined as follows. Let C F be all the vertices 

of V which are surrounded in counterclockwise sense by the circuit C* in the graph G* and let Vg, C be 
the vertices of V on the circuit C* . respectively 7^ denote the set of arcs in T going out respectively 
entering the set in G and being adjacent to one of the vertices in V^,, i.e., 

T^, = rnu.t{vi,)nu±{v^,), 
=rn^g(yi)n^g(^^.)- 

We finally define 

%{n = \*s{n + u{ri,)-Kr^.)- (1) 

The following theorem shows how the adapted reduced cost can be used to reformulate the interdiction 
problem. 

Theorem 9. Let B he a fixed budget. The following statements are equivalent. 

i) There exists an interdiction set R in G such that there exists an elementary cut [y , V \ V] in G{R) 
with u^(R\[V',V\V']) <0. 

ii) There exists a circuit C* in G* with 7|;(C*) < 0. 

^ Using the notation G instead of G* would be more consistent. To simplify notations we chose the second form. 
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Proof. Let R be an interdiction set in G and [V , V\ V] be an elementary cut in G{R) with strictly negative 
value. We partition the set R into R = VrUEr with Vr = RnV and Er = Rn E. Let [/+ = T n uj±{V') n 
u}i{VR) respectively U- = T r\uj'^{V') fl w^(Vr) the subset of arcs in T going out respectively entering into 
the set V' and being adjacent to vertices in Vr. Furthermore let W — oji {V \ R). The value of the cut 

[V, V \ V] in G{R) can be rewritten in the following way. 

^c?(fl)([y'^ V\V'])^ u{W) + uiU+) - TiU^) 

Notice that in the network G \ {RU W), there is no path from V' to vertices in V \V' . This implies 
that we can find a counterclockwise circuit C* in G* consisting only of dual arcs of arcs in Er U W and of 
arcs being adjacent to vertices in Vr. Such a circuit furthermore satisfies — V . This implies U+ = 7^ 



and U- = T~^. By a reasoning identical to the one in Section |5j we have X*g{C*) < u{W) implying finally 

Conversely let C* be a circuit in G* with 7|j(C*) < 0. We will show how to find an interdiction set R 
and a cut [V , V \ V] in G{R) with v^^^'>{[V', V \ V']) < 0. The existence of such a cut impHes directly the 
existence of an elementary cut satisfying the claim. We set V' = and Vr = V^, . As before by defining 

U+ = rC]ujt{V')r\uj±{VR) and U- = TC]ljZ{V')C]ujZ{Vr) we have U+ = T^ and U-=T~. Let C E* 

be the subset of all arcs in E* which are contained in the circuit C* , this corresponds to the set of all arcs in 
C* that are not adjacent to a vertex in Vr. Furthermore we define E^, = {e G E \ G ^^J"' minimal 

value of the cut [V, F \ V] in G{R) with respect to the interdiction set R can be reformulated as follows, 
rnin i^^^"\[V' ,V \V']) < min u^^^'^'^^'^\[V' ,V \V']) 

RCVUE:c{R)<B EnCE:c{EH)<B-c(VR) 

= 4-ciVR){Ec,) + u{U+) - ?([/_) (2) 

Note that we have furthermore \*g_-^^y^-^{E*~^) = A^(C). This comes from the facts that all arcs e* G 

C*\E'^^ satisfy A*(e*) = oo and thus have to be interdicted in the calculation of A^(C) and that c* {C*\E~^) = 
c(V/{) as already seen in Sectionjsj Further developing Equation|2]by using this relation, the definition of 7|j 
and U+ = 7^ , U- = we finally get 



min v^^"'H[V',V\V'])<j*BiC*)<0 



RCVUE:c{R)<B 

proving the claim. □ 

In the next step, we introduce a new length function x* '. E* {oo}, which is a slight adaption of A* 

and satisfies X^(C*) = 7|;(C*) for every circuit C* in G* , where x*b is the reduced length with respect to x* 
and the budget B. Such a length function allows us to solve the network security problem on G as in the 
previous sections by looking for a circuit in G* with negative reduced length with respect to x* ■ 

Let C* be a circuit in G* . For every arc e* G E* we set = A*(e*). The value of x* on the arcs in 

E^ will be defined such that for each circuit C* in G* and each vertex v E V the following equality will be 
satisfied. ^ _ _ 

Tie* n Lo^. (v)) = A*(c* n lo^, («)) + ^l(r^ n u^. (v)) ~ TiT^- n ^g. («)) (3) 



Summing the above equation over all vertices on the circuit C* gives x*(C* ni?^) = X*{C* r]E^) + u{T^J — 

1{T~), which implies the desired property Xhi^*) — 7s(^*) circuit in G* . In the following, we define 

F on the arcs EV in^ch a way that Equation|is satisfied. 

For all arcs e* G E^ adjacent to a vertex in S UT we set x*(e*) — A*(e*) — oo. This ensures that 
Equation [3] is satisfied for allv G SUT. Let v G V\{SUT), T+{v) be all arcs in T going out of v, T' be all 
arcs in T entering v, T{v) = U T~ and let be the set of all arcs in E^ adjacent to v. Furthermore, 
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let e* be an arbitrary fixed arc in . We denote by H the (geometric) graph {V*,E* U T) obtained by 
adding the arcs T to the planar graph {V*,E*). For 6^,62 G E^ , we define 7^(e^, 62) to be the set of all arcs 
in Ty that are traversed in H when going in counterclockwise sense from around w to e^. Furthermore, let 
7^+ ( 61,62) respectively 7^^(6^,62) be the arcs in 7^(6^,63), that are entering respectively going out of w. To 
simplify notations, we define 0(6^,62) = 5/(7^+ (6^,62)) — 1{T~ {el,el)). Finally we define for e* G E^ 



A*(e*) + a(e*, e*) if e* enters w 
A*(e*) — a(e*,e^) if e* leaves v . 



Let C* be a circuit in G* and v €z V a, vertex though which C* passes. We denote by ejjj respectively e*^j. 
the incoming and outgoing arc of C* with respect to the vertex v. Using the introduced notation. Equation^ 
can be rewritten as 

X*(e*n)+X*{Kut) = A*(?r„) + A*(S:ut) +a(^n,^ut) • §) 
Applying furthermore the definition of x* in Equation [3[] we finally have to prove the following proposition. 

Proposition 2. 

Sin, e„) - a(eout, e„) = a(ei„, e^^t) 

Proof. We distinguish two cases describing two possible constellations of the arcs el,e*^,e*^^. around the 
vertex v. When going in counterclockwise sense around v beginning at ejjj, these arcs are either encountered 
in the order (g'*„^e*^,t^e*) or (e?;, e*, 6*^^). 

Case 1: (ef^, e^^^, e^) 
The result follows immediately by the definition of a and the relations 

HX^ietn^Kut)) +u{T+(Sl^^,'il)) = u(7;+(e*„,e;)) 

Case 2: {e*^,el,e*^^) 
With the same reasoning as in the first case we get 

— ^ \ / — ^ / — ^ — ^ \ 

The proposition thus reduces to a(e*,e*^^) = —a(e*^^,el), which follows by observing that by construction 
of T and the fact that v ^ S UT we have 

u(r+)-r(T-) = o. 

□ 



The network security problem on G can thus be solved as in Section 6.2 by looking for a circuit C* in 
G* with negative reduced length with respect to x* ■ ^ and only if such a circuit exists, it is possible to 
diminish the value of a maximum flow in the network G by at least one unit given the budget B. It follows 
from the proof of Theorem [9] that in th^s case, an interdiction strategy can be found in the same way as in 
presented in Section [i] As the network G used in this section is only at most a constant factor larger than the 



auxiliary network used in Section 6.2 the algorithm can be implemented with a running time of 0{Bn^) or, 
by reversing the roles of budget and length, with a running time of 0{i'g'^^{G)n^). Here, too, we expect that 
the running time can be improved by exploiting more deeply the structure of the network G in the dynamic 
programming steps. 



7 Conclusions 

We proposed a planarity-preserving transformation that allows to incorporate vertex removals and vertex 
capacities in pseudo-polynomial interdiction algorithms for planar graphs. Additionally, a pseudo-polynomial 
algorithm was introduced for the problem of determining the minimum interdiction budget needed to make 
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it impossible to satisfy the demand of all sink nodes. The algorithm works on planar networks with multiple 
sources and sinks where the sum of the supplies at the source nodes equals the sum of the demands at the 
sink nodes. However this algorithm does not seem to extend easily to general network interdiction problems 
with multiple sources and sinks. Thus, it is still not known whether network flow interdiction on planar 
graphs with multiple sources and sinks is solvable in pseudo-polynomial time. We showed that the /c-densest 
subgraph problem on planar graphs can be polynomially reduced to a network flow interdiction problem on a 
planar graph with multiple sources and sinks. The algorithms presented in this paper can easily be adapted to 
the case when multiple resources are needed for removing arcs and nodes, still remaining pseudo-polynomial. 

The main purpose of the algorithms presented in this paper, was to show that various interdiction problems 
on planar graphs can be solved in pseudo-polynomial time. We expect however, that it should be possible 
to speed up the proposed algorithms by using more elaborate techniques as for example nested dissection. 
Furthermore, it would be interesting to extend the presented work to nearly planar networks as for example 
networks with a bounded crossing number or genus. 
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